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We have studied the restoration of chiral symmetry in lattice QCD at the finite 
temperature transition from hadronic matter to a quark-gluon plasma. By mea- 
suring the screening masses of fiavour singlet and non-singlet meson excitations, 
we have seen evidence that, although flavour chiral symmetry is restored at this 
transition, flavour singlet (f7(l)) axial symmetry is not. We conclude that this 
indicates that instantons continue to play an important role in the quark-gluon 
plasma phase. 



1 Introduction 

QCD thermodynamics describes hadronic/nuclear matter at finite temperature 
and/or baryon number density. This is relevant to the early universe, neutron 
stars and relativistic heavy ion collisions (RHIC). 

Hadronic matter undergoes aJxansition to a quark-gluon plasma at a tem- 
perature T — Tc ~ 100-150 MeVBcl. This transition is beheved to be a second 
order phase transition when to„, rud — > (m^ appears to be too large to make 
this transition first order). From now on we shall ignore the s quark. 

At this phase transition the SU{2) x SU{2) chiral flavour symmetry is 
restored. This means that in the plasma phase the effective (screening) masses 
of the TT and cr(/o) become equal as do the masses of the r]/r]' and 6{ao). 
The question remains as tawhether the anomalous U{l)axiai symmetry is also 
restored at this transition cla If so, all 4 screening masses would become equal 
above the transition. The rest of this talk presents lattice measurements aimed 
at resolving this question, and discussion of how this is related to topological 
charge and fermion zero modes. 

If the U{l)axiai symmetry is not restored at the chiral(deconfinement) 
transition, this means that instantons still play an important role in the plasma 
phase. 
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The functional integral for a field theory at Euclidean time in a region 
which is infinite in each of the spatial directions but has a finite extent, 1/T, 
in the time direction with periodic (antiperiodic) boundary conditions, is the 
canonical partition function for the field theory at temperature T. This allows 
us to simulate lattice QCD at finite T by simulating it on a lattice whose 
spatial dimensions are much larger than its time dimension, 

This talk presents work which is reported in detail inti, which should be 
consulted for a more complete description and for rnace complete references. 
For related work which draws similar conclusions seeQu. 

2 Meson propagators and Topological Charge at finite Tempera- 



In QCD, mesons can be created or annihilated by operators quadratic in the 
quark fields. Thus a meson propagator has the form 



For the tt, O = 75T, for the a, O — \, for the rj/rj' ^ O = 75 while for the 5, 



The connected parts of the flavour singlet and corresponding non-singlet 
meson propagators are the same. In addition to the connected part, the flavour 
singlet meson propagator has a disconnected part. These are shown schemat- 
ically in figure |I| where the gluon fields, containing the effects of closed quark 
loops have been omitted. If this disconnected contribution to the a and ri/rj' 
propagators were to vanish in the plasma phase, then the U{l)axiai symmetry 
would be restored. 

When the SU{2) x SU{2) flavour symmetry is restored. 



ture 



(1) 



O = r. 



(2) 



and 



Grj> — —Gs 
;/ symmetry were restored. 



(3) 



If the U{1) 



Gtt — —Gs 



(4) 



and 



Ga — —Gri' 



(5) 



in addition. 
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Figure 1: Connected (left) and disconnected (right) contributions to meson propagators in 
a fixed background gauge-field (not shown). 

Now let US discuss the role of instantons. These arguments are appropriate 
to the high temperature, quark-giuon plasma phase, where we can take the 
fermion mass to zero at finite volume, and we expect a gap between zero 
modes and non-zero modes at mass=0. The spectral decomposition of the 
quark propagator is 



Now let us return to the consideration of meson propagators. For configura- 
tions with no zero modes (A ^ 0), only the connected contributions survive 
the chiral (m —^ 0) limit. For configurations with one zero mode, 




S{x,y) ~ 1/m 



(7) 
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and the dominant contributions to both the connected and disconnected terms 
are ~ l/m?. Such configurations contribute with weight oc (from the 
fermion determinant), and so we get a finite contribution. The fermion deter- 
minant suppresses contributions from configurations with > 1 zero modes. We 
find Gtt = —Ga and G,,' = — G5, i.e. chiral SU{2) x SU{2) flavour symmetry 
is restored. 

Se relationship to topology is expressed by the Atiyah-Singer index the- 

m J d'^xTi-fsS = nL-nR^ Qtop = j d^xTvFF, (8) 

for a given gauge field configuration. This tells us that the topological charge 
is the difference between the number of left-handed and the number of right- 
handed zero modes. 
A corollary is 

(Q2 ) Jd^x!d^yGf{x,y) , 
Xtop = ^ = m = TO Xs (9) 

where a non-zero value of xf"^ would indicate that the 77/77' propagator differed 
from that of the tt, and that the U{\)axiai symmetry was broken. 



3 Meson propagators on the lattice. 

In lattice QCD, the SU (2) x SU{2) and U {l)axiai symmetries are explicitly bro- 
ken. We have worked with staggered fermions where 1 generator of SU{2)axiai 
remains unbroken. There are no exact zero modes, and one can only take 
TO ^ after taking the continuum limit. We must therefore work at small but 
finite m. 

We work with configurations generated in simulations on a IG'^ x 8 lattice 
in lattice QCD with 2 light (m/T — 0.05) staggered ijuark flavours. Most of 
these were provided by the HTMCGC collaborations. The chiral transition 
occurs at 5.475 < Pc = &/gl < 5.4875. 

The meson screening masses have been obtained from exponential fits to 
propagators for spatial separations. Local operators (noisy point source) were 
used for the tt, a and S. Covariant 4-link operators were used for the 7r4 and the 
r]' . t/(l) noisy estimators were used to extract the disconnected contributions 
to the a and 77' propagatorsO. 

Figure |^ shows the meson screening masses obtained from exponential fits 
to these propagators, (n.b. the flavour non-singlet scrocning masses, i.e. those 
from the connected propagators have been reported in0.) We observe that ttIo- 
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Figure 2: cr, tt and 5 screening masses as functions of /3 for full 2-flavour lattice QCD. 



and ms remain distinct until appreciably above the transition. Since the a and 
S propagators differ only by their disconnected pieces, it is these disconnected 
pieces which are responsible for this difference. The cr and tt masses are only 
distinct above the transition because of the explicit chiral symmetry breaking 
from using a finite quark mass. 

The scalar disconnected susceptibility (which would be the same as the 
pseudoscalar in the symmetric phase in the chiral limit, if we had continuum 
symmetries), shows a clear peak at the transition, and is finite well into the 
plasma phase. 
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4 The fate of zero modes on the lattice.^ 

Staggered quarks preserve one U{1) subgroup of the SU{A)axiai symmetry. 
The rest of the SU{4:) x 5'C/(4) chiral symmetry is broken by terms 0{a^). As 
in the continuum, the eigenvalues of the Dirac operator have the form: 

iX + m (10) 

and occur in complex conjugate pairs. "75" is an operator which involves 
a displacement by 4 links. It does not anticommute with the m = Dirac 
operator, and does not even have eigenvalues ±1. 
In the continuum, 

(A|75|A)-0, A^O; (11) 

and 

(0|75|0)-±1, A = 0. (12) 

For staggered quarks this is no longer true. This is fortunate since all modes 
have A 7^ 0, except on a set of configurations of measure zero. In order to 
have a sensible continuum limit, those eigenmodes with A close to zero should 
have (AI75IA) close to a constant which can be renormalized to ±1, while other 
eigenmodes should have (AI75IA) close to zero, at sufficiently small coupling. 
These near-zero eigenmodes should obey 

n+ - n_ = Qtop (13) 

where Qtop is an estimate of the topological charge of the configuration, deter- 
mined by other means. 

For each of our configurations we have determined the 16 eigenvalues A 
which are smallest in magnitude, and their corresponding eigenvectors. We 
have repeated this determination for a set of quenched configurations on the 
same size lattice. (AI75IA) was calculated for each mode. Qtop was measured 
for each configuration by the coohng method. We have plotted (AI75IA) versus 
A for the 8 smallest positive eigenvalues, A for each configuration for a repre- 
sentative (3 value for our quenched and unquenched simulations in figures ^, ^ 

For the quenched case, at /3 = = 6.2 in the high-temperature phase, it 
is clear that there is a correlation between small eigenvalues and large (AI75IA). 
In addition, we see that such eigenvalues tend to occur for configurations with 
non-zero topological charge. In the full 2-fiavour case at beta=5.55, also in 
the quark-gluon plasma phase, there is an indication of similar behaviour. 

''The discussion at the beginning of this section is condensed from Smit and Vink0 
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Figure 3: Chirality versus A at /3 = 6.2(quenched QCD). 

However, since f3 = 5.55 represents a somewhat larger lattice spacing the 
effects are far less obvious. 

We have used these low lying eigenmodes to approximate the quark prop- 
agator. Well above the phase transition this approximation gives a very good 
approximation to the disconnected part of the flavour singlet meson propaga- 
tors (sec figure H) . This supports the conclusion that the disconnected part of 
these propagators is given entirely by instantons, in the chiral limit. 

Because of the finiteness of our would be zero modes, we are forced to 
remain at finite quark mass. If we continued to zero quark mass m J d'^xTij^S 
would vanish, as would the disconnected parts of our singlet meson propaga- 
tors, leading us to the erroneous conclusion that the U{1) axial symmetry was 
restored. 
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Figure 4: Chirality versus A at /3 = 5.55(full QCD). 



5 Better Actions 

To get more definitive results would require either using much larger lattices, or 
using improved actions which better approximate the flavour chiral symmetries 
of the continuum. 

As a first try wcjJaave used the "link + staples" action suggested by the 
MILC collaboration Eil as a way of improving the flavour symmetry of the 
staggered quark action. This improved the approximation to 4-fold degeneracy 
of the eigenvalues, which would be present in the continuum. However, the 
overall improvement in the (AI75IA) versus A plot was not great. 

We are now investigating how Wilson, and Wilson + Clover Leaf actions 
perform in this regard on the configurations we have. These actions also have 
chiral symmetry breaking, but they have exact SU{2) flavour symmetry. The 
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iS = 5.55 




Figure 5: Comparison of the disconnected pseudoscalar correlators obtained from a noisy 
estimator (top curve) and from the truncated spectral decomposition of the quark propagator 
(bottom curve) at /3 = 5.55 {Nf = 2,m = 0.00625). 



way which they fail to obey the Atiyah-Singer index theorem is different from 
the way staggered fermions fail. 

Finally, we are investigating how well domain-wall fermions satisfy the 
Atiyah-Singer index theorem and looking at their zero modes on these con- 
figurations. Domain-wall fermions are 5-dimensional Wilson fermions with 
boundary conditions in the 5th dimension which ca.use light quark states to 
reside close to the boundaries in the 5-th dimensionE3. Their advantage is that 
in the limit of infinite 5-th dimension, chiral SU{2) x SU{2) flavour symmetry 
becomes exact, and the approach to this chiral symmetry appears to be ex- 
ponential in N^. In this limit, the index theorem should become exact. Our 
preliminary calculations on a few configurations show promise. 
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Other improved actions, including ones with gauge improvement should 
be tried. 

6 Discussion and Conclusions 

At the phase transition from hadronic matter to a quark-gluon plasma in QCD 
with 2 zero mass quarks, the chiral SU{2) x SU (2) chiral flavour symmetry is 
restored. We have found evidence that the anomalous U{l)axiai symmetry is 
not restored, but remains broken in the plasma phase. The explicit breaking of 
chiral symmetry by the staggered lattice fermions is sufficiently large at these 
lattice spacings to somewhat obscure these results, as does the explicit chiral 
symmetry breaking due to our use of finite quark masses. We are prevented 
from decreasing the quark masses to zero by this chiral symmetry breaking 
which renders would-be zero modes non-zero. Since our analytic arguments 
required a fixed spatial volume, we need to repeat our simulations on several 
spatial volumes to study finite size effects, to see that our conclusions survive 
in the infinite volume limit. 

One can improve on our results by setting the 4Qfop eigenvalues of smallest 
magnitude to zero by hand. This gives us a non-zero m / d'^xTrj^S and non- 
zero disconnected contributions to our singlet meson propagators. However, 
the chiral flavour symmetry breaking associated with staggered fermions mani- 
fests itself in ways other than the zero mode shift, and these are not addressed 
by this procedure. To address these requires us to study more systematic 
improvements to finite temperature simulations and measurements involving 
light quarks. 

If the U{l)axiai symmetry remains broken in the high temperature phase, 
this means that instantons continue to play an important role in this phase. 
In particular they contribute the disconnected part of the propagator which is 
responsible for the mass difference between isosinglet and isotriplet scalar and 
psciidoscalar mesonic excitations of the plasma. In addition, they make con- 
tributions to the connected meson propagators. This occurs because, despite 
the fact that the fermion determinant vanishes as in the chiral limit for the 
Qtop = il sector, the connected and disconnected contributions to the meson 
propagators for gauge configurations in this sector diverge as giving a 

net finite contribution. Thus, while = m„ and m^/ = m^, ttIt^ ^ m^/ 
and ttIct ^ mg. It will be interesting to see if the meson propagators in the 
Qtop — sector yield screening masses close to 27rT, the screening mass for 
a free quark-antiquark pair. If so, it would indicate that it is the instantons 
which prevent the quark-gluon plasma from being described as a gas of free 
quarks and gluons. 
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It is important that the chiral flavour symmetry is SU{2) x SU{2). For 
5'C/(3) X SU{3) chiral flavour symmetry, instantons are further suppressed in 
the plasma phase, and the U{l)axiai symmetry is restored at the transition, 
which now becomes first order. 

At high temperatures, the disconnected contributions to the flavour singlet 
meson propagators are well approximated by the contributions of a few low 
lying modes of the quark propagator. 

We need improved actions — especially fermion actions — to get more 
definitive results, unless we have the computing power to work with the much 
larger lattices needed to enable us to decrease the lattice spacing enough to sig- 
nificantly reduce the lattice breaking of chiral flavour symmetry. Domain- wall 
fermions appear to be a good bet. It is of interest to know if such improve- 
ments also give one better critical indices at the phase transition. Others have 
suggested that the "perfect action" improvements are the most promisingEj. 

Even with the ultimate improved action, currently used simulation meth- 
ods will fair poorly when used to calculate meson propagators in the chiral 
limit of the plasma phase. This is because, as mentioned above, conflgurations 
which occur with probability only O(m^), each contribute a term 0(1/™^) 
to the meson propagators. Hence, in the chiral limit, important, and in fact 
inflnite contributions come from conflgurations which never occur. We are 
studying algorithm modiflcations to better accommodate this limit. 
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